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Abstract. We introduce a ^-deformation of the genus one s?2 Knizhnik-Zamolodchikov- 
Bernard heat equation. We show that this equation for the dependence on the moduli 
of elliptic curves is compatible with the qKZB equations, which give the dependence 
on the marked points. 



1. Introduction 

The Knizhnik-Zamolodchikov-Bernard equations are a system of differential equa- 
tions arising in conformal field theory on Riemann surfaces. For each g,n G Z> , a 
simple complex Lie algebra g , n highest weight g -modules Vi, and a complex parameter 
k, we have such a system of equations. In the case case of genus g — 1, they have the 
form 

(1) Kd z .v = - J2 hi 3) d K v + - z h \)V>%. 

The unknown function v takes values in the zero weight space V[0] = r\ x< z^Ker(x) of the 
tensor product V — V\ ® ■ ■ • (g) V n with respect to a Cartan subalgebra f} of q . It depends 
on variables zi, . . . , z n G C and A = ^ X u h u G I) , where (h u ) is an orthonormal basis of 
P) , with respect to a fixed invariant bilinear form. The notation x^> for x G End(Vj) or 
x G g means 1® ■ ■ • ®x® ■ • ■ ®1. Similarly x™' denotes the action on the ith and jth 
factor of x G End(V r i ® Vj). 

The "r-matrix" r G g <8> g obeys 

r{z, A) + r(-z, A) (2il) = 0, [r(z, A), h ® 1 + 1 <8> h] = 0, V/i G fj , 
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with (^2 i Xi®yi)( 2lS) = J2i Ui® x ii an d is a solution of the classical dynamical Yang-Baxter 
equation jFWf {r^^ = r(zi - z 2 , A) <g> 1 G £/g 03 etc.) 

V V V 

_[ r (1,2)^(1,3)1 _ ^(1,2)^(2,3)1 _ ^(1,3)^(2,3)1 = q 

As a consequence, the KZB equations ([I]) are compatible, meaning that if the equations 
are written as VjV = 0, then the differential operators Vj commute with each other. 
The solutions of the classical dynamical Yang-Baxter equation arising in conformal field 
theory are parametrized by the modulus r in the upper half plane and can be expressed 
in terms of theta functions, see IIKWL IF VII . 



A difference version of this story was proposed in [[FJ: Suppose that for an Abelian 
complex Lie algebra I) we have h -modules Vi, i = l,...,n with a weight decompo- 
sition Vi = ©|te{j*Vj[//] into finite dimensional weight spaces Vi[/i]. Then we say that 
meromorphic functions Rij(z,X) of z G C and A G t) * with values in Endf, (Vi <g> Vj), 
(1 < i 7^ j < n) form a system of dynamical R-matrices if they obey the (quantum) 
dynamical Yang-Baxter equation 

Rij(zi -z 2 ,X- 2r ] h (3 Y 2) R l k(zi ~ zs, X) {13) R jk (z 2 -z 3 ,X- 2 V h^Y 23) 
= R jk (z 2 - z 3 , X) {23) R ik (zi -z 3 ,X- 2 V hW)( 13) Rij(zi ~ *2, ^) (12) , 
in End(Vi <8> Vj ® V&) for alH < j < k and are "unitary": 

Rij(z, X)Rji(-z, XY 21) = ldvi®v r 

We adopt a standard notation: for instance, R(z, X^r/h^)^ acts on a tensor Vi®v 2 ^>v 3 
as R(z, X — 2?7/i 3 ) <g> Id if v 3 has weight /x 3 . 

The deformation parameter ("Planck's constant") is here rj. If we have a family of 
dynamical .R-matrices depending on rj such that Ry = Idy^v- + 2>Wij + 0{rf) as 77 — > 0, 
we recover the classical dynamical Yang-Baxter equation and the unitarity condition for 
Tij, i. e. the properties that r obeys, viewed as an element of End(\^ (g) Vj). 

If we have a system of dynamical .R-matrices Rij we can then construct a compatible 
system of difference equations, the qKZB equations for a function v(zi, . . . , z n , A) G V[0]. 
They are a dynamical version of the I. Frenkel-Reshetikhin qKZ equations |[FR]| , and 



their semiclassical limit are the KZB equations. Their construction is reviewed in |2Tl 
below. 

The main examples of solutions of the classical and of the quantum dynamical Yang- 
Baxter equations are associated with elliptic curves. In the quantum case, they can be 
viewed as intertwining operators between tensor products of representations of elliptic 
quantum groups taken in different orders |[FV2|| . In the rank one case (one-dimensional 



f) ) explicit expressions for R matrices -Ra,m depending on two "highest weights" A, M G C 
are known. They are associated to pairs of evaluation Verma modules [ |KV2 | for the 



elliptic quantum group E TjV (sl 2 ) and were computed using the functional realization of 
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these modules [[FTVl]] . If n highest weights Ai, . . . , A n G C are given, then Rij = Ra^Aj 
form a system of dynamical R matrices as described above. 

Hypergeometric solutions of the corresponding qKZB equations were introduced and 
studied in |[FTV1|| , |[FTV2|| , ||FV4|| . See also [TJ where similar equations are studied and 



solved. Special cases of these equations appear in the statistical mechanics of RSOS 
models. Form factors and correlation functions in the infinite volume limit are con- 
jectured to obey qKZB equations. In these cases explicit formulae were proposed by 
Lukyanov and Pugai |[LPf |. 

The subject of this paper is a deformation of the KZB heat equation: in the classical 
case, additionally to the KZB equations above, that are associated to the variation of the 
marked points on the elliptic curve, one also has an equation associated to the variation 
of the modulus r of the elliptic curve. The function v also depends on r and one has an 
additional equation, compatible with the KZB equations, the KZB heat equation 



AniKd T v = A\v + - s(z, A, r) 



{ij) v. 



l :J 



for some s G g <8> • Here A a denotes the Laplacian of f} corresponding to the invariant 
bilinear form. For example, if n = 1 then this equation reduces to 

AnindrV = (A A - ^ p(«W> r)e a e_ a )v, 

where e a are properly normalized root vectors and p is the Weierstrass function. 

In this paper we propose a discrete version of the KZB heat equation in the rank 
one case. The heat operator is an integral operator, whose kernel is given in terms of 
hypergeometric integral solutions of the qKZB equations of |[b'TV2|| . 



In Sect. we review the qKZB equations and their hypergeometric solutions. Then we 
introduce the elliptic Shapovalov form, which is an ingredient in the integral operator, 
and the qKZB heat equation in Sect. [| We prove that it is compatible with the qKZB 
equations, discuss its properties and show in Sect. f|, in an illustrative example, that its 
semiclassical limit coincides with the KZB heat equation. Finally, in Sect. |5] we study 
the special case where the step of the difference equation is a negative integer multiple of 
the deformation parameter. In this case, the semiclassical limit gives the KZB equations 
with positive integer k, the situation arising in conformal field theory. In this case the 
the KZB equations are defined on sections of the finite dimensional vector bundle of 
conformal blocks, of which we describe a difference analogue in simple cases. 

It is very likely that the hypergeometric solutions of the qKZB equations are also 
solutions of the qKZB heat equation. However, we were able to prove this only in 
the case where the sum of the highest weights is two. In this case the hypergeometric 
integrals are one-dimensional. 

In a sequel to this paper, we show that integral operators of the kind introduced in this 
paper can also be used to describe the transformation properties of hypergeometric so- 
lutions under the modular group. In fact it turns out that the hypergeometric solutions, 
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at least if A; = 2 obey remarkable identities under transformations of the modulus r 
and the step p by SL(3, Z) acting on CP 2 with affine coordinates r, p. These identities 
give both the solutions and the monodromy of the solutions. The whole picture results 
in an non-Abelian version of the properties of the elliptic gamma functions, which is a 
generalized Jacobi modular form for SL(3, Z) x Z 3 in the sense of [ FYS] . 



Acknowledgment. We thank R. Ferretti for explanations on Gauss sums. 



2. Hypergeometric solutions of the qKZB equations 

2.1. The qKZB equations. Fix A = (Ai, . . . , A n ) G W 1 such that m = J27=i A i/ 2 is a 
nonnegative integer, and a complex number rj. Unless stated otherwise, we will assume 
that these parameters are generic. 

Let t and p be generic complex numbers in the upper half plane. 

Let V^. be the vector space with basis e ,ei,... equipped with the action of an 
operator h given by he^ = (Aj — 2k)ek- We view as a representation of the Abelian 
Lie algebra f) = Ch. 

To these data is associated a system of dynamical i?-matrices and thus a system 
of qKZB difference equations. The i?-matrices R\.\ k {z, A, r) [ |FV2| | of E TtV (sl<i) are 
endomorphisms of VV ® V& k . Let = V\ 1 ® • • • <g) V& n . The qKZB equations are 
equations for a meromorphic function v(z, A) of z G C n and A G C taking its values in 
the zero weight subspace Vjj[0] = Ker(^™ =1 of (this subspace is nontrivial since 
^2 Aj/2 is assumed to be a nonnegative integer). It will be more convenient to view v as 
a function v(z) taking values in the space JF(V^[0]) of meromorphic functions of A G C 
with values in Vr[0]. Let 5j, j — 1, . . . , n be the standard basis of C n . Then the qKZB 
equations have the form 

v(z + p5i) = Ki(z,T,p)v(z), i = 1, . . . ,n. 

The qKZB operators Ki(z, r,p) act on the space ^"(V^[0]) and are given by 

Kj(z, T,p) = Rjj-^Zj - Zj-i +p, T) ■ ■ ■ Rj,l(Zj - Z!+p, t) 

TjR jtn (zj — z n , t) ■ ■ • Rjj + i(zj — z j+ i, t). 
The operators Rj t k{z,r) are defined by the formula 

fe-i 

R jtk (z,T)v(\) = R AjAk (z,\-2r) h«\r)v(\), 

and (Tjv)(X) = v(X — 2r)n) if h^'v(X) = fxv(\) and is extended by linearity to JF(V^[0]). 
The qKZB system of difference equations is compatible, i.e., we have 

(2) Kj(z+pSi,T,p)Ki(z,r,p) = K t (z + p5 j ,T,p)K j (z 1 T 1 p), 
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for all j, I, as a consequence of the dynamical Yang-Baxter equations satisfied by the 
i?-matrices. We also consider the "mirror" qKZB operators 



Kj & Pi T ) = R j,j+l( Z j ~ Z J+l + T > P) ■ ■ ■ R j,n( Z 3 -Zn + T, p) 



TjR^Azj - z u p) ■ ••R)j j _ 1 (z j - Zj- X ,p), 



with 



Rl k (z,p)v(\) = R AjAk (z,\-2r) ^ ^,p)v(X), 

The corresponding system of qKZB equations 

v(z + r5j) = Kj(z,p,r)v(z), j = 1, . . .ra, 

is also compatible. In fact, if we write x y = (x n , . . . , Xi) for any x = (x±, . . . , x n ) G C™ 
and let P : V% — > V^ v be the linear map sending v± <S> ■ • • ® v n to v n ® • • • ® v i, then we 
have, adding the dependence on A in the notation, 

KVftp, r; A) = P^K^&tP, r; A v )P. 



2.2. Hypergeometric solutions. In |[FTV2|| we constructed a "universal hypergeo- 
metric function", which is a projective solution of the qKZB equations: it is a func- 
tion u(z, A, fMjTjp), defined for generic values of the parameters 77, A, taking values in 
Vjj[0] Cg) Vjj[0] and obeying the equations 

w (£ + <^p, r, p) = Kj (z, t, p) ® Dj u(z,T,p), 
(3) u(z + 5jT,T,p) = D) ® KJ(z,p,t)u(z,t,p), 

u(z + 5j,T,p) = u(z,T,p). 

Here we view w as taking values in the space of functions of A and fi with values in 
V^[0] (g> V^[0]. Aj acts on the variable A and Kj on the variable fx. The operators 
Dj, act by multiplication by diagonal matrices Dj(fx), DJ(X), respectively. For our 
purpose, the most convenient description of these matrices is in terms of the function 

a(X) = exp(— niX 2 /Arj). 

We have, for j = 1, . . . , n, 

D .U) = °(f-MhM + ■■ ■ + eT¥jEc;Ar£Lj Al) 

AN a(ji-2ri(h® + '-- + h<ri)) 

D v (x) = «(A - 2n{hW + ■■■ + hW)) M ^ Al „y^ , Al) 
iK) a(X - 2r/(M 1 ) + • • • + /#))) 

These operators are diagonal in the basis of Vj[0] formed by tensor products ej = 
® • • • ® ej n of basis vectors of the Va*. so that XX^fc — 2i&) = 0. 

From w one can construct projective solutions (eigenfunctions of the corresponding 
difference operators) by taking coefficients of the basis vectors ej. If D i {fx)ej = dij(fx)ei, 



dij(fj,) G C, and u = J2 uI ® e i then for any fixed / and fi, the function v(z,X) = 
u 1 (z, A, n,r,p) obeys v(z + p5i) = d i j( y fj,)K i (z,T,p)v(z). It follows that 

n 

v{z,X) =]Jd i>I (ji)-^ p v(lX), 

i=i 

is a true solution to the qKZB equations. The parameters I and fi determine the 
multipliers, as is easily seen from the explicit expression for u below: 

(4) v(z + 5i, A) = d itI (/i)~ l/p v(z, A), v(z,X + l)=e ^ d(z, A). 

The second system of equations in (^) gives the monodromy of these solutions, see 
|FTV2J. 

The explicit expression for u is given by the following formulas. 

(5) u(z,X,fi,r,p) = e~^r / J^O^ Afe {U - Zk,r,p) Y\_^-2 V (ti - tj,r,p) 

i,k i<j 

}]uJi(t, z, A, r)uj(t, z,fi,p)dti ■ ■ -dt m ei <g> ej. 
i,j 

The phase function Q has the product formula 

'I — e 2Tri(z-a+jT+kp)\M _ e 2-Ki(-z-a+(j+l)T+(k+l)p) \ 



j,k=Q ^ 



<>■-■ '■!■>■■- g| e 27rt(z+o+ir+fcp)Vl _ e 2m{-z+a+(j+l)T+(k+l)p)\ 



It is symmetric under exchanging r and p and obeys the functional equation 



(6) n a (z + p, r,p) = e 2m J{z ± °' T) Q a (z, r,p). 

The weight functions uji are given by 

u + \ \ Tf Q(ti-tj,T) ^ A-p-r ^j- 6{U - z k + r]A k ,r) 

u ° w <*■■■■• n ^-, +2ff „) ,£ n n n 

n-j-r 0(tj - tj + 2r], t) -A- -j-r 9(X+tj-z k - i]A k + 2r]i k - 2rj YllZl ( A f - 2i f ) , r) 
H 0(t f . r ) 1111 0(L - z fc - nAfc, r) 

k<uei k ,jeii Kl 31 ' k=ijei k Vj fe ' fe ' y 

The summation is over all n-tuples I\, . . . , I n of disjoint subsets of {1, ... , m} such that 
I k has i k elements, 1 < k < n. The theta function is here the first Jacobi theta function 

The "mirror" weight functions are related to the weight functions with the reversed 
order of factor. Indicating the dependence on the highest weights explicitly, we have 
ujj(t, z, ji,p, A) = dj/v(t, z^, p,,p, A v ), where, as above, (x\, . . . , x„) v = (x n , . . . , x\). 
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The integral over t\, . . . ,t m of the 1-periodic integrand in fl5|) is defined by analytic 
continuation from a region of the space of parameter where the rule J^Aj = 2m does 
not hold: one starts from the region Im(r/) < 0, Im(r), lm(p), Im(^Aj) > for which 
the integral is over the torus (M/Z) m and defines the integral in general by analytic 
continuation. 



2.3. Remark. In |[FTV2|| we used only weight functions and no mirror weight functions. 



Then the qKZB equations (|3|) only involve qKZB operator and no mirror qKZB oper- 
ators. The choices of this paper make the qKZB heat equation more transparent. The 
proof that u obeys the relations ([3]) is the same as the proof of Theorem 31 in JFTV2J. 
Note however that the conventions in the definitions of Dj are different there. 

3. The qKZB heat equation 

In this section we define a q-analogue of the KZB heat equation, prove that it is 
compatible with the other qKZB equations and show that the differential KZB heat 
equation arises in the semiclassical limit in the simplest non-trivial case. 

The qKZB heat equation is an integral equation. The integration kernel is a contrac- 
tion with the fundamental hypergeometric solution. The contraction is defined using 
the elliptic Shapovalov form. 

3.1. The elliptic Shapovalov form. For j = l,...,n, /x, r 6 C, Imr > 0, let 

Q Kj (/!, r) : ®VL — > C be the bilinear form on with matrix elements Q Kj (/i, r) (e^® 

ei) = 5 k ^ j {fi,T), 

(7) nV s - ( r ) V TT 9{2r]{A i + 1 ~ °' T)6{27]1 > T) 

U ^*^ T) \9(2 V ,t)J l\e(n + 2r l (A j + l-k-l),T)9(ii-2r l l,T)- 

Out of these bilinear forms we define a bilinear form Q(/j,, t) : (g) — > C on the 
tensor product 

n-l 

Qfa T ) = Q Al (fji, r) Q A2 {fi + 2 V hV , r) ( 2 - n+2 ) ■ ■ ■ Q An {fi + 2i] ^ h® , r) ^ , 

and a bilinear form on the space of functions of A with values in Vjj[0]: if / and g are 
holomorphic functions from C to Vjf[0], we set (if the integral converges) 

QAf ®g) = J Q(^, T)f(n) <8> g(-n)a(fji) dfx. 

The integration is on the path t i— > 2vjt + e, — oo < t < oo. This bilinear form is called 
the elliptic Shapovalov form. 

The main property of the elliptic Shapovalov form is that the i?-matrix is in a cer- 



tain sense symmetric with respect to it, see Lemma 3.9. In particular it is a sort of 



contravariant form for the action of the elliptic quantum group, see eq. ( |T2"D below. 
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3.2. Notation. To write the following formulae in the most transparent form we will 
use the following notational conventions. Let for k £ Z>i and a complex vector 
space V, Fkiy) denote the space of meromorphic functions of k complex variables 
with values in V® k . For example u(z,r,p) £ JF 2 (V^[0]). We also set J r o{V) = C 
and F{V) = Fi{V). If / £ ^(V) and g £ F k (V), we define f ® g e F j+k by 
(/ ® #)(Ai, . . . , Aj+fc) = /(Ai, . . . , Xj) ® g(X j+1 , . . . , A i+Jfe ). If ^ £ End(F(V)) are dif- 
ference operators with meromorphic coefficients we write A\ ® • • • ® A r £ End(J>(V)) 
to denote the composition of the operators Ai, each acting on the zth variable and 
the zth factor. We also use this notation if one of the Aj is an integral operator 
Q : D C F 2 (V) -> C of the form / ^ / $(//)/(//, Q(p) £ (V ® V)*, de- 
fined on some subset -D. Then Ax <g> • • • ® A r is defined on some subset of J- r {V) and 
maps to JF r _ 2 (y). 

3.3. The qKZB heat equation. Let T(£, r, p) be the integral operator on V^[0]-valued 
functions of one complex variable A 

(8) T(z, t, p)v = (a® Q T+P )u(z, t,t + p) ® v, 

where a is the operator of multiplication by the function a(X). More explicitly, if {ej} is 
a basis of V^[0] consisting of tensor products of basis vectors and u = J2i j u i,J e i ® ej, 
v = J2 v i e i> T )( e i ® e j) = Qi(^ r )^,J ; we have 
(9) 

T(z,T,p)v(X) =^ ^a(X)^2 j u l! j(z,\,fi,T,T + p)Qj(fi,T + p)vj(-fi)a(fi)df?j e 7 . 

Theorem 3.1. The equations 

v(z + pdj,T) = Kj(z,T,p)v(z,T), j = l,...,n, 

(10) 

v(z,t) = T{z,T,p)v(z,T +p), 

are compatible, i.e., we have, in addition to 

T(z + p5 j ,T,p)K j (z,r + p,p) = Kj(z,T,p)T(z,T,p). 



The proof of this theorem is given in |3.6| . 

A similar statement holds for the qKZB difference operators KJ: they obey the com- 
patibility identities 

KY(z + t6i,p,t)K?(z,p,t) = K^(z + T5j,p,T)Kj(z,p,T), 

for all j, I, and the operator 

T y (z,p,r)v = (Q T+P ® a)v ®u(z,t + p,p), 
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obeys 

T v (z + Td j} p } T)Kj(z,p+T } T) = KJ(z,p,t)T v (z,p,t), 

implying the compatibility of the mirror qKZB equations 

v(z + r8j,p) = Kj(z,p,r)v(z,p), j = l,...,n, 
v(z,p) = T s/ (z,p,r)v(z,p + t). 

Corollary 3.2. Let U(z, X, p,r,p) be the function 

U(z,r,p) = a <g> Q T+P ® a(u(z, r,r + p) <g> u(z, r + p,p)). 
Then U obeys the system of equations jj^). 

Conjecture 3.3. Let U(z,r,p) be the function defined in Corollary \3.2j . Then 

U(z,r,p) = Cu(z,T,p) } 

for some constant C. 

This conjecture is proved in the case where J^Aj = 2 (with C = —e imr) / (27T\Amt)). 
This proof will be published elsewhere, ||FV3|| . 

Assuming the conjecture correct, we can obtain solutions to the full system ( |10D as in 
|373| : for arbitrary I — (h, . . . ,i n ) G Z> with ^(A^ — 2ik) = and \i G C, let 

2 n 

v(z,\,t) = YldiM'^u^z, X,n,T,p). 

i=l 

Then the function v(z,t) : A i— >■ v(z,X,r) is a solutions of (pT0|) . It also obeys (|j) and 
f (£, A, r + 1) = e i ^t>(i', A, r). 

Remark. The Shapovalov pairing Q T contains an integration J dp, which we chose to 
be the integral on the path t i— > 2r]t + e. This choice makes the integral convergent if 
Im(?7) < for a large class of functions, thanks to the strong decay at infinity of the 
Gaussian function a(p) on this path. This class contains in particular our hypergeomet- 
ric solutions. 

It should be however emphasized that the only properties of J dp that are needed for 
this construction are that it be a linear form on functions of p invariant under translations 
by 2r) times weights of vectors in V\. n and that it be well defined on a suitable class of 
functions. 

In particular, if the highest weights are rational with greatest common denominator 
d, we may replace the integral over p by the sum over the set {A + 2rjk/d, k G Z}, so 
that the heat equation may be viewed as a difference equation of infinite order. 
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3.4. The case of integer highest weights. If A G Z> , let 5a be the [) -submodule 
of Va generated by e^+i, ga+2, ■ ■ ■ ■ The A + 1-dimensional quotient V\/ S\ is denoted 
by L\. The classes e~o, . . . , of e , . . . , ca build a basis of L\. The space La carries 
a one-dimensional family of representations of the elliptic quantum group E TV (sl2), see 
FV2 |. 

For integer highest weights A iy Aj, the .R-matrix R\ u \.(z, \,t) preserves S\ i ® Va 7 



and V\ t (g> S'a |[F V2| , |[FTV1]| . Therefore it induces an endomorphism, still denoted 



Rhi,Kj{. z i A, r), of La^ <S> L^.. If Ai, . . . , A n G Z> we then have a system of dynamical 
i?-matrices and thus a system of qKZB equations defined on L^[0] = (L^®- ■ -®La„)[0]. 

A universal hypergeometric function A, /i, r, 77) taking values in the tensor product 
of finite dimensional modules L^[0] ® L^[0] and obeying (|3|) was found in |[MV|| . It is 



defined by A, /x, r, 77) = jUij(z, A, /i, r,p, r])e/ ® ej, where uj t j are the analytic 
continuation of the components of the universal hypergeometric function for [0] <S> [0] 
which are shown to exist for these values of /, J. 

Then we can introduce a heat operator T(z,r,p) acting on functions with values in 
Lfi[0] by the same formula (g). 

Theorem 3.4. Suppose that A\, . . . ,A n are non-negative integers. Then the equations 
v(z + p5j,r) = Kj(z,T,p)v(z,r), j = l,...,n, 

(11) 

v(z,r) = T(z,t,p)v(z,t + p), 
for a function v taking values in L^[0] are compatible, i.e., we have, in addition to (f|j ; 
f(z + p5 j ,T,p)K j (z,T + p,p) = Kj(z,T,p)f(z,T,p). 

The proof of this Theorem is contained in |3. 7| below. 

3.5. Rational 77. A particularly interesting case is the case of integer highest weights 
and rational 77. Let us for instance assume that 2Nr] = 1 for some positive integer N 
and suppose that the highest weights A 1; . . . , A n are positive integers. 

If N is large enough, then the qKZB operators may still be defined. Indeed we have: 

Lemma 3.5. Let Ai, A2 be positive integers, and N be a large enough integer. Then the 
matrix elements of the R-matrix Ra 1 ,a 2 { z i A, r; 77) G End(LA x <E> L\ 2 ) with respect to the 
basis {ei (8> <%} are regular functions of 77 at 2r) — 1/N for fixed generic values of z,X,t. 

Proof: The R matrix Rai,a 2 { z i ~ z 2, A, r; 77), for generic r], is uniquely determined up 
to normalization by an intertwining condition for tensor products of .EV^s^-modules 
LiA.(zi), see ||FV2|| . The E TyV (sl2) module La (z) may be realized for integer A as a 



symmetric tensor product of two-dimensional modules [F V2|| , so that the matrix elements 
of RAi,Aj, for a basis consisting of symmetrized tensor products of basis vectors, can be 
expressed as polynomials in the matrix elements of Rn- The latter matrix elements 



11 



are known explicitly and are regular as functions of rj. For i = 1,...,A, the basis 
vector gj of L\ is proportional to the symmetrized tensor products of basis vectors 
of the two dimensional modules. The proportionality constant is an elliptic factorial 
E[j=i Qfivfy/Qfiv) which is regular and non-zero at 2r]N = 1 as long as iV > A. 

Thus if N > max(Ai, . . . , A n ), the matrix elements of the R- matrix are regular at 
2r]N = 1. □ 

Fix some generic complex number e. Then we may consider the qKZB equations as 
equations for functions v(z, A), where the dynamical parameter A runs over the finite set 
{k/N + e\k G Z/2A^Z}. Indeed, the coefficients of the qKZB operators are 1-periodic 
functions of A and the shifts of A in the difference operators Tj are integer multiples of 
2rj = 1/N. The shift by the generic number e ensures that on this finite set no poles of 
the qKZB operators are encountered. Thus, we get: 

Proposition 3.6. Suppose that N = (2r])~ 1 and A\, . . . ,A n are positive integers, with 
N large enough. Fix a generic complex number e. Let -F/v(e) be the space of functions 
f : iZ -> Vj[0] so that f(X+2) = /(A). Then the qKZB operators K^z, r,p), K?(z,p, r) 
are well-defined endomorphisms of F^{e). 

In this situation we thus have a truly holonomic system, i.e., a compatible system 
of difference equations for functions taking values in a finite dimensional vector space 
Fjv(e). In order to define the heat equation we have to worry about the fact that the 
universal hypergeometric function u is not defined for all values of the parameters. Recall 
that u(z, A, p,,T,p) is also a meromorphic function of r\. Let us say that r\ is a regular 
point for u if u is regular at this point for all A, fi and all generic z, r,p. 



Theorem 3.7. Let rj, A = (A 1; . . . , A n ), e be as in Proposition \3.o\ and assume that rj 
is a regular point for u. Then the heat operator 

2N-1 2 

T N (z,T,p)v(X) = e — (1 ® Q(ii k ,T + p))u(z, A, /i h ,r,T + p) ®«(-/i t )e * , 

with fik = —e + k/N, maps Fjy(e) to itself. Moreover, the equations for v(z, r) G F^(e) 
v(z + p5j,T) = Kj(z,T,p)v(z,T), j = l,...,n, 

v(z,t) = f N (z,T,p)v(z,T + p), 

are compatible, i.e., we have, in addition to ^), 

f N (z + p5 j ,T,p)K j (z,T + p,p) = Kj&T^fn&TiP), 

on F N U). 



The proof of this Theorem is contained in p.7| below. 

The description of the set of regular points for u will be studied elsewhere. Here we 
only remark that in the case n — 1, Ai = 2, the point rj = 1/2N is a regular point for 
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all N > 3, as can easily be checked since u is given by a one-dimensional integral. In 
this case, Conjecture |3]^ holds, see ||F V3|| , namely, we have 



u(z, A, v, T,p) = C N e 2" 



27V- 1 



X 



(1 ® Q(nk, T + p)® l)(u(z, A, /i fc , r, r + p) <g> i/, r + p, p))e * 



fc=0 



for all A G e + 4Z, z/ G -^Z. The constant is C 



AT 



; e 2Tri/N 
S(N) ' 



with the Gauss sum 



27V-1 



S(N) = = (1-*)V^V. 

£i=0 



3.6. Proof of Theorem |3.1| . The proof is based on some identities involving R- 

matrices, Q and Dj. As above, we set a(A) = exp(— m\ 2 /Arf) 



Lemma 3.8. For any A, M 7 



a 



(A - 2r]h( 2 )) 



R AM (z + r, A, r) = e-^^iWz, A, r 



Proof: One way to prove this lemma is to use the functional realization (see [b'TVl]): 
The matrix -Ra,m relate two bases of the same space of functions. The basis elements 
are products of ratios of theta functions and have therefore well-behaved transformation 
properties under shifts of z by r. The computation is straightforward and will not be 
reproduced here. □ 



Lemma 3.9. Let A, M G C and v , w G V A ® Vm- Then 

(Q A (/i + 2 V h^ 2 \ r) ® Q M (/i, r) u, i2 AiM (z, r) w) = 

, t)v, w). 

Proof: We first prove a version of this identity for L-operators. Let L((, A) G End(C 2 <8> 
V\(z)) be the L-operator of the evaluation Verma module V\(z). We claim that, for any 
v u v 2 G C 2 ® V A (s), 

(QH/x + 2r//i (2) , r) ® QV r) Vlj L(C, -//, r) u 2 > = 
(12) = (Q\fi, t) <g> Q A (fi + 2 V h m , t)L(C, \i + 2 V (h^ + /i (2) ), t)v u v 2 ). 

We have Ql(fj,,r) = 1 and Q\([i,t) = 6{2rj)6(n — 2r])~ 1 9(fi)~ 1 . Define the matrix el- 
ements of L by L((,fi)ej ® v = J2k=oi e k ® L*j((,fi)v. Then the claim is equivalent 
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to 



Qk(v, T )( e k, L° ((, -fi)e k ) 

= Q A (/i + 2 V , r)(L°(C, fi + 2 V (A - 2k + l))e k , e k ), 
Qk(P» r)(e h , L°(C, -/i)e fe _i) 

= Qj(^ ^gt^/i - 2rj, r){L l (C, fi + 2^(A - 2k + l))e fc , e^}, 
+ 2r/(A - 2*), r)Q A (/i, r)(e k , L (C, -fj)e k+1 ) 

= Qt +1 (fM + 2r), r)(L?(C, fx + 2 V (A -2k- l))e k , e k+1 ), 
Q\Qi + 2 V (A - 2k), r)Q^, r)(e k , L\{C -fJt)e k ) 

= QI(jm, t)Q£(h - 2rj, t)(L\(C, fi + 2 V (A -2k- l))e fc , e k ). 



These identities follow immediately from the explicit expressions given in [|FV2|| for the 
operators L k (called a, b, c, d in |[FV2 



We now extend this result to the general case. We use the intertwining property of 
the i?-matrix: let C\,Cm be the L-operators of V\(zx), Vm(z 2 ), respectively. ThenQ 
Ra,m( z ~l ~ z 2i £ End(VA(zi) ® Vm(z 2 )) is uniquely determined up to a factor by the 
relation 

The operators C and £ (giving the action of the quantum group on the tensor product 
by using the coproduct and the opposite coproduct, respectively) are defined by 

£(C,a*) = £ A (C,/i-2^ 3 ))( 12 )/:M(C^) (13) , 

C'(C,fi) = C M ((,ix-2r,hW)WC A ((,^ 12 \ 

The .R-matrix normalized by the condition Ra ; m(zi — z 2 , /x)eo <8> eo = eo <%> e$. 

In particular, if v = w = eo ® eo, the claim of the lemma is correct for trivial reasons. 
We prove the general case by induction: let us suppose that the lemma is proved for 
v, w of weight A + M — 2j, j = 0, . . . , k — 1, k > 1. Now it is known, see [ |FV2| |, that, for 
generic parameters, the weight space V\(zi) ®Vm(z 2 )[A + M — 2k] is spanned by vectors 
of the form £°(C, A) x (or C\{C, A) x), ( G C, x of weight A + M — 2(k — 1), and any fixed 
generic A. Indeed, if these vectors did not span the weight space, they would be part of 
a proper submodule, contradicting the irreducibility of the tensor product. 

By iterating (|X2|) , we obtain 



(Q\fi + 2 V (h® + h®)) ® Q A (/x + 2 v h^) <g> Q M (/iK C{(, 
= (Q\fi) ® Q A (/x + 2t?(/i« + ^)) <g> Q M (/x + 2^«) 
£'(C, /i + 2r ? (/ i ( 1 ) + + /i (3) )K, w 2 ). 



-//)v 2 ) 



1 Here we omit the argument t to shorten the notation 
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In particular, if v\ — eo ® v, v<i — e\ <S> w, one has 

(Q A (y + 2 V h®) ® Q M (/i) v, £?(C, -/i) to) 
= Q{(n)(Q A (v + 2t7(-1 + h^)) <g> g M (^ - 2^)£'J(C, /i + 2tj(l + + /i (2) )) u, w). 

We turn to the proof of the induction step. Suppose that v, w have weight A + M — 2k, 
and write w = C'1((,n)x. Let us set z = z\ — z 2 . Then 

(Q A (fi + 2r)h®) Q M (fi) v, R AM {z, -y) w) 
= (Q A (y + 2 V hW) ® Q M (y) v, R AM (z, -rfC'liC, -//) x) 
= (Q A (fi + 2 V h®) <g> Q M (/i) v, CUC, -v)Ra,m(z, -V + 2 v) x) 
= Q\(fi)(Q A (fi + 2 V (-1 + &<*>)) ® g M (/i - 2 V )£'l((, n + 277(1 + fc« + /»<»))) v, 

= Q\(y)(Q A (y - 2rj) ® Q M ( / u + 2r/(-l + fcW)) 

#A,m(*, /i + 2r7(-l + h w + /i (2) ))/:'J(C, /i + 277(1 + A + M - 2A;)) v, x). 

In the last step, we used the induction hypothesis. The calculation continues by com- 
muting R with £', and then by bringing C to the right. This last part is similar to 
the above calculation read backwards, and will not be reproduced in detail. One finally 
obtains, as desired, 

(Q A (y)®Q M (y + 27 1 h (1) )R A 

This completes the induction step and thus the proof of the lemma. □ 



Lemma 3.10. 



a(D])- l K j (z,T,p + T) = K j (z,T,p)ae-™ r,A i ( Z w *i) 



a Dj Kj (z, p,r + p) = Kj (z, p,r)ae 



Proof: This is a straightforward consequence of the definition of the difference operators 



Kj , Kj and of Lemma |3J|. □ 



Lemma 3.11. Let f , g be holomorphic functions from C to V\. 

Q T+p (f, Kj{z, r + p,p)g) = Q r+p {Dj l K](z + p5 3 , r + p, r)f, g) e ^^w^ 
Q T+p (Kj(z,r + p,r)f,g) = Q^f^D^Ktf+rS^T 



p,p)g)e 



mriAj(J2ijij A;) 



Proof: The proof of the first identity is given by using Lemma |3.9| to bring the R- 
matrices in Kj to the left, the translation invariance of the integral to bring Tj to the 
left, and |3.10| to commute the resulting Kj (z + p5j, r + p, —p) with a. The proof of the 
second identity is similar. □ 
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We can now complete the proof of Theorem |3~TL Let Cj = e mvAi ^ 1 ^ and v a 
function from C to V^[0]. 

T(z+p5 j ,T,p)K j (z,T + p,p)v = 
= (a® Q T+P )u(z + p5j,r,p + r) <g> Kj(z,t + p,p) v 
= Cj(a® Q T+P )(1 <g> Dj 1 Kj(z + p5 j ,T + p,r)u(z + p5j,r,p + t)) ®v 
= C j {a®Q T+p ){{D y j Y l ®Dj l u{z + (p + r)^,r,r + p)) <g> v 
= Cj{a ® Qr+^dD^^Kjiz, t, t + p) ® 1 r, r + p)) ® v 
= Kj(z,T,p)T(z,T,p)v. 



3.7. Proof of Theorem |3.4| and Theorem |3.7| . Theorem |3.4j can be proven in the 
same way as Theorem FO. 



There is however an apparent difficulty: the heat operator involves the Shapovalov 
form which contains a sum over all components of u, including those that a priori do not 
have a limit for integer highest weights. The solution is provided by Theorem 2 of ||MV||: 



let us say that I = . . . , i n ) G Z> is admissible for A if i a < A a for all a = 1, . . . , n. 
Then (a special case of) Theorem 2 states that the components ujj(z, A, /x, r,p) such 
that / or J is admissible, are regular functions of the highest weights at A for generic 
values of the other variables. Moreover, Q A (fi,r) vanishes if k > A, cf. (|7|), so that 
the sum appearing in the compatibility condition is effectively restricted to admissible 
indices. 



Theorem |3.7| is proven in the same way as Theorem |3J] and |3.4j . In fact the only 
property of the integral over \x that is used in the proof is the translation invariance. So 
the same proof gives the compatibility relation in this case provided the function of /i^ on 
the right-hand side is periodic in k with period 2A^. Now u(z, X, /i, r,p) is exp(— ■niNXfi) 
times a 2-periodic function of A and fi. So the exponential factors combine into the 
expression 



iTTN(\ + tl k ) 2 



If A G e + jjZi, this expression is periodic in k with period 2A^. The same argument 
shows that Tn(z,t,p)v(\) is 2-periodic in A for A G e + -^Z. □ 



4. Semiclassical limit 

We consider here the semiclassical limit of our quantum heat equation in the simplest 
non-trivial case and show that we do recover the KZB heat equation in this limit. The 
case we consider is n = 1, with Ai = 2. The qKZB equations for the dependence of z\ 
are trivial in this case and we can assume that z\ = 0. The zero weight space is one- 
dimensional, and we identify it with C using the basis e±. Suppose that v v (X, r) is a fam- 
ily of solutions of the qKZB equations with parameters r,p = —2kt],t,T], parametrized 
by rj around zero. Assume that v v has an asymptotic expansion v v (X, r) = v (X, r)+0{rj) 
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Figure 1. The integration cycle 7. The points ±2r) are the singularities 
of the integrand 

at 77 = 0. We want to find the equations satisfied by Vq. For this we expand the qKZB 
heat equation 



(13) 
v v (X,r) 



around 77 = 0, setting p = —2kt] and keeping r, k, A fixed. The dependence of 77 of the 
constant in front of the integral was chosen in such a way that the semiclassical limit 
exists. The integration path is t *■ p, — rjt (t G R). The hypergeometric solution u is 
independent of z in this case and is given by the formula: 

M \ [ r, u ,6(\ + t,T) 6(p + t,T + p) Jm 

u [\,fi,T,p,r] ) = e 2" il 2ri {t,T,T + p)— — -— — ■ — -dt. 

J<y 0(t-2ri,T)0(t-2ri,T + p) 



The integration cycle 7 is depicted in Fig. |l|. 



Theorem 4.1. Suppose that v v (X, r) is a family of solutions of ftTBj) with an asymptotic 
expansion v v (X, r) = Vq(X, t) + r)Vi(X, r) + ■ • • , then v(X, r) = v (X, r)/9(X, r) obeys the 
KZB heat equation 



dv 

2niK-— 



d 2 v 



Or dX 2 



2p(X,r)v + c{r)v, 



for some c(r) independent of X. 
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Proof: The integral on the right-hand side of ([13]) has the form 



/Air) 



r/oo 



rjoo 



This integral has the asymptotic expansion as rj —>■ 



g(\,\,0)+r) 



1 d 2 



in dfi 2 



d 

g(X,/i, 0) + — 
fi=x 9 V 



g(X,X,r)) )+0(r, 2 ). 



r)=0 



To compute the various terms of this expression, we first notice that the integration 
cycle in u is pinched by the singularities as r) —>■ 0. The integral defining u can then be 
expressed as a divergent (as r) — > 0) part given by 2ni times the residue at t = 2r) plus 
the integral on a cycle 7 which stays away from the singularities. 
To compute the residue we introduce fl by 



n, 2 r,(t,T,T + p) 



I _ e 2ni(t-2 V ) 

2iri 

= (t - 2r))-^—— [ n 2rl (2r), r,r + p) + 0((t - 2 V ) 2 ). 

As 7) — > 0, fi 2 ^(2?7, t,t — 2nr)) is regular and converges to 1. 
We then have 

2?rc n (n , v 0(A + 277, r + p)6{4r), t + p)v v {(i, r + p) 
-^—jil 2r) (Zr], t,t + p) 



g(x,^,v) 



1 



2ixi 



x ■ 



9'(0,r)9(fx + 2r],T+p) 

fi 2r? (t, r, r + p)— — - — — —-dt 

9(t — 2r),T) 9{t - 2t],t +p) 

9{Ar),T + p)9'{0,r + p) 



9(p + 2r),T + p)9{)i -2 V ,t + pf^ T + P) - 
From these formulae we can compute the various terms: 

g(X,X,0) = v (X,t), 

d 2 „ ... „„ _o fv (X,r) 



Finally 



d_ 

dr) 



9 2 /i 
g(X,X,r)) 



g(\,fi,0) = 9(,\,T)di 



0(X,t) 



d(r) v (X, t) + t) 



r)=0 



d_ 

dr) 



v v (X,t) 



r)=0 



2 

Til 



9{t + \,r)6(t- X,t) 
9(t,r) 2 9(X,r) 2 



dt v (X, t). 
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Here Ci(r) is some scalar function independent of A. Using the identity 
9(t + X,r)9(t- X,r) 1 



9(t,r) 2 9(X,r) 2 ~~ 9'(0,r) 2 
we see that the right-hand side of (O) is 



(p(A,r) - p(t,r)) 



v v (X,t) 

r]=0 

Id 2 d 2 \ «o(A,r) 2 . 



w (A,r) +r/— 

for some function c(r) independent of A. 

Since the first two terms also appear on the left-hand side, the proof is complete. □ 

5. CONFORMAL BLOCKS 

In this section, we introduce, in the simplest case of one marked point, a difference 
analogue the vector bundle of conformal blocks. We begin by reviewing the differential 
case. The vector bundle of conformal blocks is, in this vector bundle on the 

moduli space A4i,i of genus one curves with one marked point. The projectivization of 
this vector bundle carries a connection given by the KZB differential operator. We then 
give a difference analogue of this vector bundle. It has a (discrete) connection, which is 
now given by the qKZB heat operator T. 

5.1. The differential case. Let g be a simple complex Lie algebra with Cartan subal- 
gebra I) and root space decomposition g = h©© ae A0a- Let the non-degenerate invariant 
bilinear form ( , ) on g ~ g* be normalized so that the highest root 9 obeys (9, 9) = 2. 
Let k be an integer larger than or equal to the dual Coxeter number h y of g, and A e f)* 
be a dominant integral weight, so that (9, A) < k — h v . Denote by L\ the irreducible 
g -module of highest weight A. 

To these data one associates a holomorphic vector bundle of conformal blocks on 
the moduli space Ai\,i of genus one complex curves with one marked point |TUY|| . Its 



projectivization carries a canonical flat connection. The fiber over a point may be defined 
as a space of coinvariants for the Lie algebra of g- valued rational functions on the curve 
whose poles are at the marked point, acting on the irreducible affine Kac-Moody Lie 
algebra module of highest weight A and level k — h y . 

An explicit description |[FW|1 of this bundle, which for our purposes can be taken as 



a definition, may be obtained by viewing M.\^ as the quotient of the upper half plane 
H + by SL(2, Z). We may then regard the vector bundle E k> a of conformal blocks as an 
SL(2, Z)-equivariant vector bundle over H + . Let La[0] = {v E L\\t)v = 0} be the zero 
weight space space of La- It carries a natural linear action of the Weyl group W of g. 
The fiber of E Kj a over r G H + is then to the space of holomorphic maps v : I) — > La[0] 
such that 
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(i) v(X + qi + q 2 r) = exp(— 7ri«(g 2 , <?2)t - 2mn{q 2 , X))v (A), for all A G f) and gi, g 2 in 
the coroot lattice Q v . 

(ii) v(w ■ X) = e(w)w ■ v(X) for all w G W , where e : W i— > {±1} is the homomorphism 
sending reflections to —1. 

(iii) For all roots a, a; G a , and integers I > 0, r, s, the map t> obeys the vanishing 
condition 



as a(A) — > r + sr. 

The action of SL(2,Z) on the base may be lifted to an action on the bundle: let 



obeying the cocycle condition ipghij) = i>g{h ■ r)^(r). Denote by n : H + — > -Mi,i 
the canonical projection. Local holomorphic sections of the vector bundle of conformal 
blocks on an open set U C are then the same as holomorphic sections v of E KtTn on 
7r _1 (?7) so that v (g ■ r) = ip g (T)~ l v (r). In other words, they are holomorphic functions 
v(X, t) on C x obeying (i)-(iii) for each fixed r and such that 



The projectivization of this vector bundle carries a holomorphic connection, and hor- 
izontal sections may be constructed by an elliptic version of hypergeometric integrals 



We describe here the connection in the case of si (2, C). If g = si (2, C) and A = ma, 
rn — 0,1, ... , then L\[0] is one dimensional. Let us chose a basis of La[0] and identify t) ~ 
f)* with C via the basis a/2. Then E k ^(t) = E K ^m[j) consists of holomorphic functions 
v(X) on the complex plane so that (i) v(X + 2r + 2st) = exp(—27iin(s 2 r + sX))v(X), 
(ii) v(—X) = (— l) m+1 f(A), (iii) v is divisible by 6(X,T) m+1 in the ring of holomorphic 
functions. 

If k > 2m + 2, we have £ Kj2m (r) = ^(A, r) m+1 K _ 2m _ 2 (r) w/ , where Q k {t) w is the 
k + 1-dimensional space of holomorphic even functions obeying (i). Otherwise £? K)2m (r) 
is trivial. 

It follows that 



x l v (X) =0((«(A) 



r - sr) l+1 ) 





FV1I. 
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The connection on E K}2m is denned by its covariant derivative T(U, E K)2m ) — * r(C/, E K ^m)® 
on local holomorphic sections: 

Vv(X, r) = (d T - -L (d 2 x - m{m + l)p(A, r)) - 71^)-%^) J v(X, r) dr. 
Here p is the Weierstrass elliptic function with periods 1 and r and 

oo 

v ( T ) = e niT/12 Y[{l-e 2nijT ) 
i=l 

is the Dedekind //-function.^ In spite of the poles of the p function, this connection is 
well-defined on E K ^ 2m as can be seen by noticing that the poles cancel in the expression 
of the induced connection #~ m ~ 1 o V o Q m+1 on Sj^l 2m _ 2 . The fact that V preserves (i) 
and (ii) is easily checked. 

The connection V is SL(2, Z)-equivariant, in the following sense: if U C H + is an 
SL(2, Z)-invariant open set, and g G SL(2, Z), we have the pull-back g* : T(U, E K ^m) —* 
T(U, E Kt 2 m ), sending a section v(t) to ^ £ ,(r) _1 t>((7-r). We may extend g* to T(U, E K ^ m )® 
by tensoring with the pull-back of differential forms. Then j*oV = Voj*. 
Therefore the connection is well-defined on the vector bundle of conformal blocks on 
Mi tl . 

Example. If m = 0, V is essentially the differential operator of the heat equation. The 
theta functions 

rGZ+j/2re 

form a basis of K (r) for fixed r, and obey the heat equation 2nind T 6j yK = d\6j, K . 
Moreover, we have 0j tK (— X, r) = 9-j tK (X,r). It follows that the functions 

(14) Vj(X, r) = ^(r)- 1 {e j+1>K (X, r) - 0-,_ M (A, r)), j = 0, 1, . . . , « - 2, 

form a basis of the space of horizontal sections. 



See [ FV1 | for the case of arbitrary m. 



5.2. The difference case. Let us turn to the difference case (for si (2, C)). We describe 
a difference analogue of E^ 2m , a holomorphic vector bundle E^ 2m ,ri on H + which is 
preserved by the qKZB heat operator. We fix a generic rj in the lower half plane. 
Guided by the semiclassical analysis of Sect. H, we suppose that —p/2r] = k is an integer 
> 2 and consider the qKZB heat operator (|8|) for n = 1, z% — 0, = 2m. 

We start with the somewhat trivial but instructive case m — 0, and write T k>0 (t) = 
T(z — 0, r,p — —2r)K). Here the qKZB heat operator is 

T hJ o(t)v(\)= [ e-^ x+ ^v(-^dfx. 



2 The connection, being on the projectivization, is really defined up to adding a multiple of the 
identity. We have chosen it here so that it defines a connection on the vector bundle over Mis 
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The integral is over the path 1 1— > 2r/t, —oo<t<oo. 

We define E Kj2 m=o,r] = E Kt0 to be the holomorphic vector bundle of odd theta functions, 
as in the differential case: the fiber over r G H + is E KOrj (r) — {f G @ k (t) | /(—A) = 
-/(A)} 

Theorem 5.1. Let k > 2 and suppose that Imr] < 0, Imr > 0. Taen T Kj o(r) maps 
E K , ,r,(r - 2t]k) to E KfitV (r). 

This theorem is based on the identity 

6 j)K (\, t) = -!= [ e-^+^e^i-ft, r - 2r)K) dp, j G Z/2kZ, 

V 4z7 7 J2r]R 

which gives the action of T Kj0 (r) on the basis 9j — 9_j, j = 1, . . . , k — 1, of 6 K (r — 2t\k). 

Let us now turn to the case of general m. To compare with the classical limit we 
consider the qKZB operator for the quotient v of the dependent function by YlT=i ^(^ + 
2?7j, r), i.e., we set 

T«, ot (t) = m (r) _1 o T(z = 0,r,p = r - 27/k) o m (r - 2t]k), 

where m (r) is the operator of multiplication by the function A i— > njli ^(^ + 2nj, r). 

Example. If m = 1, the qKZB operator for v is v ^ T K i(r)v is 



^ , k,i( t )' u (A) = ct(^) / V(\, /j,,t,t — 2T]K)a(fi)v(—[i)diJ,, 

J2TpL 

with kernel 



TviXj.1 j 

V(X, /i,r, a) = ce Q 2r) (t,r,a 



6(t - 2r], t)6(X + 2% r)6(t - 2rj, cr)6(n + 2rj, a) 
for some c = c(r, a) independent of A, /i. The integration cycle is depicted in Fig. [I]. 

Let E Kj2 m,r]( T ) be the space of holomorphic functions so that 

(i) v (A + 2r + 2sr) = exp(47u?7m(m + l)s — 2niK,(s 2 T + sX))v (A), 

(ii) ^(-A) = (-ir+ i n7 =1 fg±^^(A), 

(iii) f is divisible by rij=o ^(^ — ^Tjj, t) in the ring of holomorphic functions. 

Alternatively (and more simply), E K ^m,r l {T') is the space of functions of the form 
Y\^ =Q 0{X - 2ryj,T) tp(X), with ip G 6 K _ 2m -2(T) w/ . In particular, E Ky2 m, v {r) has the 
same dimension as the space E K ^m{^) appearing in the differential case. Let E K ^m,r] — 
2m,T]( T )- It is naturally a holomorphic vector bundle over H + . 

Theorem 5.2. Let m, k G Z> , k > 2m + 2 and suppose that mi?] < 0, Imr > 0. Then 
(r ) maps E R 

,2m, 77 (j" 

- 2r)K) to E Kt 2m, v (r). 
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Proof: This theorem is a corollary of the results of [|FV4 |. We give here the proof in the 
simplest case m — 1. The proof of the general case is similar. Let v G E Kj 2 >v (t — ^V K )-> 
and set v = T Kj i(t)v. Properties (i), (ii) for v can be checked straightforwardly, by using 
the identities 

6{\ + 2,t) = 6(\,t), 9(X + 2t,t) =e- 4wi{x+r) 6{\,T), 9(-X, r) = -9(X, r), 

obeyed by 9 and translating the integration variable in the integral over fi. The latter 
involves moving the integration contour, which presents no problem as the vanishing 
condition (iii) for v guarantees that the integrand has no poles. Let us check that v is 
holomorphic and obeys (iii). As the zeros of 9(X, r) are simple and on the lattice Z + rZ, 
v is regular except possibly for simple poles at —2r] + Z + rZ. 

We claim that v vanishes at A = r + sr and at A = 2r] + r + sr for all r, s G Z. Then 
(ii) implies that v is regular at the points — 2^ + Z + rZ (and thus everywhere), and that 
v is divisible by 9(X, t)9(X — 2rj, r). 

Since v obeys (i), it is sufficient to prove the claim for r,s£ {0, 1} 
It follows from (ii) that v(0) = and, in conjunction with (i), also v(r + sr) = 0, 
r, s G {0, ±1}. For example, we have 

6(T + 2r,,r) e -^9(-r + 2 V , r) _ 8 

v(—t) = —. rv It) = —, : v(t) = — e 'v(t). 

K ' 9{t-2t],t) y ' 9(t-2t],t) k j k j 

On the other hand, (i) implies v(— r) = e~ 8mv v(T), so v(t) = 0. 

Let us check that v(2i]) vanishes. By using the functional equation (^) for Q 2 r), we 
obtain 

J 9(t-2r],T)9(Ar],T)9(t-2r],a)9(/j, + 2r],a) 
ce"^ J n 2v (t + a,T,a)^-. _ + ^ rdt 



ce"" 



9(Arj,T)9(t-2r],a)9(fM + 2rj, a) 
9(p, + t — a, a) 
9{Ar), T)9(t -2t]- a, a)9(fi + 2//. a ) 

M *' T ' a) 9(4r],T)9(t-2 V ,o-)9(fi + 2 V ,o:) dt 



= 9^) TeSx ^ V{X >^ a) - 

In this calculation the change of variable t i— > t — a was used. For this our choice of 
t-integration contour is essential, since it implies that one does not encounter poles when 
one deformes it back to the original position. For general m this identity is part III of 
Theorem 26 in [ KV4 |. Thus 

9{Ar),T) 
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But it follows from (ii) that 

so v(2rj) = 0. The same argument may be applied to 2r\ + r + sr with r,s£ {0, ±1} (or 
even for general r, s) . We have 



V (2rj + r + sr, r, a) = ) ' \ e 2mS(T ies x= . 2ri+r+ST V(X, //, r, a) 



9(Arj,r 
This implies that 

On the other hand, using (ii) and (i) we obtain the same equation but with the opposite 
sign, so that both sides vanish. Thus v(2ri + r + st) = and v is regular at the potential 
singularities A = — 2r] + r + sr, r, s G Z. □ 

A more direct reformulation of this theorem is the following. 

Corollary 5.3. Let m, k G 1i>o, k > 2m + 2 and suppose that Imrj < 0, Imr > 0. Let, 
for t G C m , 

m( ' ' j ~ m - t, + 2 V , r) jl 9( tj - 2 V m, r) " 

Introduce the integral kernel 

e V N no(X,(jL,T,p)9(fji,p) 



Lj=-m 

where 



/m 
Y[^2 V m{ti,r,p) Yl ^-2 V (ti-tj,r,p) 



i=l l<i<j<m 

xu m (t, A,r) u m (t,[i,p) dti ■ ■ -dt T 



The integration is over a torus as in 
Then the integral operator 



M(r)0(A) = / M(A,//,r,r-2t7K)0(-/i)d/i 

J2rjR 

maps Q K - 2m -2(r - 2t]k) w to Q K _ 2 m-2(r) w . 
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5.3. Remark. A section v of E Ki 2m,n is called projectively horizontal if it obeys the 
qKZB equation T K>jn (r)w(r — 2t)k) = C{r) v{r) up to a scalar factor C(r). For m = 
projectively horizontal sections are given by odd theta functions as in the differential 
case, see (|b|). In a sequel |[FV3|| to this paper, we show that for m = 1 (and conjec- 



turally for higher m as well), projectively horizontal sections are again given by elliptic 
hypergeometric integrals. 

5.4. Remark. The compatibility of the difference operator T Kjm (r) with the SL(2,Z) 
action can be better understood in terms of a discrete connection on a space with an 
SL(3, Z)-action. This will be discussed in [FV3]. 
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